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SUMMARY 


A large space structure with members missing was investigated using a 
finite-element analysis. The particular structural configuration was the tetra- 
hedral truss, with attention restricted to one of its planar faces. Initially 
the fini te-element model of a complete face was verified by comparing it with 
known results for some basic loadings. Then an analysis was made of the struc- 
ture with members near the center removed. Some calculations were made on the 
influence of the mesh size of a structure containing a hexagonal hole, and an 
analysis was also made of a structure with a rigid hexagonal insert. In 
general, load- concentration effects in these trusses were significantly lower 
than classical stress-concentration effects in an infinitely wide isotropic 
plate with a circular hole or circular rigid inclusion, although larger effects 
were obtained when a hole extended over several rings of elements. 


INTRODUCTION 

Large, low-mass structures are expected to play an important role in future 
space missions. These missions include activities such as solar power collec- 
tion, navigation. Earth resource surveillance, and communications. A discussion 
of potential space activities is presented in reference 1, and a survey of 
expected requirements of large space structures is presented in reference 2. 
Although specific missions are not yet defined, some form of large, open truss 
structure is expected to be a strong candidate for providing a stiff skeletal 
frame upon which functions can be conducted in orbit. One such structure is 
a tetrahedral truss, or octetruss, sketched in figure 1. The tetrahedral 
truss, constructed of identical column members, has psuedo-isotropic elastic 
properties. Some design considerations for large space structures, such as the 
tetrahedral truss, are reported in references 3 and 4. 

The present report examines the ability of this type of truss to function 
safely if one or more of the column members are damaged and unable to carry a 
load, a potential design problem not considered in reference 3 or 4. Such 
damage could result, for instance, from impact by meteoroids or space debris, 
assembly accidents, or docking accidents. In addition, effects that exceed- 
ingly stiff equipment might have when attached to the truss are investigated. 

A fini te-element model of a tetrahedral truss with members missing from its 
isogrid planar faces or with a rigid insert simulating the stiff equipment was 
subjected to in-plane loadings. Load concentration induced as a consequence 
of the missing members or the rigid insert was investigated only in a two- 
dimensional framework. Since redundancy is added by the core members for the 
three-dimensional truss, the two-dimensional results tend to overestimate loads 
and are conservative for design of the three-dimensional truss. 



The analytical tool used to obtain the results was the NASTRAN 1 computer 
program (ref. 5). Initially, attention was focused on a reference model of the 
two-dimensional structure with no missing members. To verify the model, compar- 
isons were made with closed- form solutions of reference 3 for basic loading 
conditions of extension in each of two directions and shear. A representative 
portion of a large mesh with selected members removed was subsequently analyzed 
under the same extensional and shear loadings to obtain load-concentration fac- 
tors. Load-concentration factors exhibited by a structure with missing members 
that formed a hexagonal hole are compared with classical stress-concentration 
factors in an infinitely wide isotropic plate with a circular hole. The effect 
of increasing the number of elements in the mesh with a constant-size hexagonal 
hole was determined. Load-concentration factors exhibited by a structure with 
a rigid hexagonal insert were also computed and are compared with stress- 
concentration factors in an isotropic plate with a rigid circular inclusion. 
Finally, the effect of increasing the number of members removed from the truss 
was examined. 


SYMBOLS 

H length of member (see fig. 2) 

n number of members along an edge (see fig. 2) 

N x ,Ny direct loading per unit length in x- and y-directions (see fig. 3) 

N X y shear loading per unit length (see fig. 3) 

Pqo,P 5 oO,P_ 60 ° loads in members inclined 0°, 60°, and -60° with respect 

to x-direction 

r^ "radius" of hexagonal hole or radius of circle that circumscribes 

hexagon 

x,y coordinate directions defined in figure 2 


REFERENCE STRUCTURAL MDDEL AND LOADINGS 
Structural Model 

As indicated previously, this study considers a planar face of a tetrahedral 
truss. The model of a planar face is an n x n mesh shown in figure 2, where 
n is always an even number. The model is a two-dimensional assemblage of iden- 
tical members (rod elements in the finite- element ccmputer program) with pin 
joints assumed at the nodes. Constraints were chosen to restrain rigid-body 
motion without inducing additional loads into the members. For the full model, 
motion in both x- and y-directions was restrained at the node halfway between 


1 NASTRAN: Registered trademark of the National Aeronautics and Space 

Administration. 
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the top and bottom of the left edge, and motion in the y-direction only was 
restrained at the corresponding node on the right edge. Therefore, rigid- body 
rotation was restrained. Equivalent boundary conditions, except for uniform 
translation in the x-direction, were imposed on quarter- symmetry models of the 
structure. 


Loadings and Load Distributions for Reference Model 

Internal member loads due to anticipated loadings in space are very low. 
Thermal loading and inertial loadings are likely to be determining factors for 
internal loads, but depend a great deal on equipment which is mounted to the 
truss. To aid in the understanding of the behavior of these large space struc- 
tures without unduly complicating the present study, some basic loading condi- 
tions were considered. 

The loading conditions, depicted in figure 3, are an essentially uniform 
loading in the x-direction N x , a shear loading N X y, a loading in the 
y-direction Ny, and a loading in the y-direction with some side loadings, 

Ny corrected. The loading N x is not quite uniform since loads of equal 
magnitude were applied at all joints along the jagged edges, even at the top 
and bottom joints. The uniform loading with side loadings (fig. 3(d)) simu- 
lates the combined loading on a portion of a mesh infinite in both directions. 
This loading is referred to as Ny corrected. Discussion of this loading 
condition is presented subsequently. 

Loading in x-direction .- For the reference model, individual member loads 
for the loading N x , as taken from reference 3 and reproduced by the present 
infinite-element analysis, are 


p 0° {3 

^ ~ ~2 


( 1 ) 


p 60° _ 
N/” ' ° 


( 2 ) 


p -60° 
= 0 


(3) 


where P 0 o, P 60 o» and p -60° are the loads in members parallel, inclined 60°, 
and inclined -60° to the x-direction, respectively. The x-direction is shown 
in figure 2. The chosen loading was applied by means of the system of concen- 
trated loads shown in figure 3(a). Each of these loads had a magnitude of 
{'JJ/ 2)Nj£ , including those at the top and bottom joints. 
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Shear loading .- For the uniform loading N X y, the individual member loads, 
as taken from reference 3 and reproduced numerically by this study, are 


P 0 o 

5x? 


(4) 


p 60° 

N x / 


(5) 


p -60° _ 
NjjyJl 


(6) 


The uniform loading was applied by means of the concentrated loads shown in 
figure 3(b). Here loads of magnitude N X y£ were applied to joints on the top 

and bottom rows and loads of magnitude \[3 N X y^ were applied to the protruding 
joints on the jagged sides. At all corner joints, loads of half the magnitude 
of loads on the corresponding sides were applied. 

Loading in y-direct ion. - For the uniform loading Ny on a doubly infinite 
mesh the individual member loads from reference 3 are 


p 0° \/3 

N y & 6 

p 60° {3 

NyS- 3 

p -60° \[3 

Nyfi- 3 

To reproduce these results numerically, special attention must be paid to the 
jagged edges. In figure 3(c), no loads on the jagged edges were considered, 
and in figure 3(d), side loads were considered. In both cases, loads of magni- 
tude Ny& were applied at the joints on the top and bottom rows except at the 
corners where the magnitude was (1/2) N y &. 


(7) 


( 8 ) 


(9) 
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Results for a quarter panel of the configuration of figure 3(c) f which did 
not include any loads on the jagged edges, are shown in figure 4 for a 16 x 16 
full mesh. The results show edge effects; that is, member loads are signifi- 
cantly different from those of equations (7), (8), and (9) near the jagged 
edges, but tend toward these values near the center. The maximum member load 
is 27 percent larger than that predicted by equation (8). This edge effect 
occurs even when more members are added to the mesh. Note that in figure 4, 
as in all figures in this paper that display nondime ns ional member loads. Nil 
was chosen to be 1 000 for convenience. The plot was produced by making a modi- 
fication to the ccmputer program described in reference 6. 

For the loading, Ny corrected, edge loadings of magnitude (VJ/6)Ny£,, 
which represent the contribution of those members that would exist in a doubly 
infinite mesh, were applied at the joints along the jagged edges, as indicated 
in figure 3(d). Member loads for the configuration of figure 3(d) agree with 
those of equations (7), (8), and (9). 


RESULTS AND DISCUSSION 

The member loads obtained in this study were nondimensionalized so that 
they are independent of Young's modulus, cross-sectional area, and length of 
the members, providing that the assumption of pinned joints remains valid. 
Models with members removed from the structure were analyzed. Two types of 
single members were removed: (1 ) a horizontal member near the center and 
(2) a diagonal member near the center. These two situations cover all possi- 
bilities of a single member being removed at the center of the truss. A model 
having an entire six-member unit about a central joint removed and a model 
having a center joint made rigid with the six-member unit made inextensible 
were also studied. Finally, additional results on models with more extensive 
hexagonal holes were obtained. 


Models With One Member Missing 

The effects of removing a horizontal member in a 1 6 x 1 6 mesh for the load-' 
ings, N x , Ny, and Ny corrected are shown in figures 5(a), 5(b), and 5(c), 
respectively. No results are shown for a shear loading because there were no 
changes in member loads from those predicted by equations (4) , (5) , and (6) . 

As can be seen from figure 5(a), the member loads are almost symmetric about 
the center of the removed member. Notice that the effect of the removed member 
dissipates pretty rapidly, as indicated by the small values of load on the diag- 
onals more than three member lengths from the removed member. A comparison of 
the two types of loading in the y-direction (figs. 5(b) and 5(c)) reveals essen- 
tially the same distribution of load except within three member lengths of the 
sides where the previously mentioned edge effects occur (compare figs. 4 
and 5 (b) ) . 

Next, the 16x16 mesh with a diagonal member removed was subjected to the 
four basic loading conditions. Member loads for the loadings N xy and Ny are 


h 
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shown in figures 6(a) and 6(b), respectively. No results are shown for the 
loading N x because there are no changes in member loads from those predicted 
by equations (1), (2), and (3). No results are presented for the Ny corrected 

because of the similarity (except for edge effects) to those displayed in 
figure 6(b). 

Maximum nondimens ional member loads in the vicinity of the removed member 
(either horizontal or diagonal) are listed in table I. Presented also in this 
table is the load- concentration factor which is defined here to be the ratio 
of the maximum load in a member near the damaged area to the load in that same 
member for the undamaged condition. With a horizontal member removed, maximum 
member loads are 26 percent and 17 percent higher than loads in the reference 
model for extensional loadings in x- and y-directions, respectively. With a 
diagonal member removed, maximum member loads are 35 percent and 26 percent 
higher than loads in the reference model for shear loading and loading in the 
y-direction, respectively. 


Models With Six Members Missing 

An entire six-member unit about the center node of a 16 x 16 mesh was 
removed to form a hexagonal hole, and the remaining structure was analyzed. The 
resulting maximum nondimens ional member loads in the structure are presented in 
table I. The maximum loads increased by 30 percent and 47 percent for exten- 
sional loadings in the x- and y-directions, respectively, and by 56 percent for 
the shear loading. The graphic displays of all nondimens ional member loads for 
three of the loading conditions are presented in figure 7. Because of the sym- 
metry of the structure and either symmetry or antisymmetry of the loadings, only 
the upper right quadrant of the structure is shown. 

Since the unbroken mesh exhibits isotropic behavior (see ref. 3 or 7), 
it is interesting to compare the present results with solutions of an infi- 
nitely wide isotropic plate with a circular hole at its center. As indicated in 
table I for extensional loadings, the classical plate stress-concentration fac- 
tors would be 3.0 for extension in the x- and y-directions, and 4.0 for shear 
(see ref. 8) . Thus the present load- concentration factors for the mesh are 
significantly lower than might be presumed. However, the two situations are 
not quite comparable since the location of maximum load for the mesh does not 
always correspond to the location of maximum stress concentration in the iso- 
tropic plate. For example, in the case of the shear loading, the maximum load 
in the mesh is in a diagonal member that joins a node on the boundary of the 
hole (fig. 7(b)), whereas classical stress concentration for shear loading is 
located on the edge of the hole. 

A procedure for designing isogrid structures (the form of mesh considered 
here) with a cover plate with a hexagonal hole has been given in reference 7. 

If this procedure is attempted when no cover plate exists for large space struc- 
tures, the results would be overly conservative since the procedure is based on 
classical stress-concentration effects for extension of an isotropic plate. 

Additional studies were made on the model with a six-member unit removed 
to see whether any significant changes would occur as the size of the structure 
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increased. In these studies, the length of the members, and hence the size of 
the hole, were constant with the overall dimensions of the structure growing 
as the number of elements in the mesh along an edge was increased. Table II 
presents the results for various mesh sizes which indicate a slight decrease 
in maximum nondimens ional loads and thus a slight decrease in load-concentration 
factors as the number of members in the mesh was increased. Plots of these data 
as a function of the reciprocal of the square of the mesh number are presented 
in figure 8 and are extrapolated to give some idea of convergence. 


Model With Rigid Insert 

A 16 x 16 mesh structure with a rigid insert was examined. The rigid 
insert was formed by rigidly connecting the center node to the adjacent six 
nodes and permitting no relative rotation between members at the center joint. 
Such a situation might simulate a piece of equipment being connected to the 
center and six adjacent joints. The solution was obtained with the MPC (multi- 
point constraint) feature in the finite- element analysis. The graphic displays 
of nondimens ional member loads are presented in figure 9. Again only one quad- 
rant of the structure is shown. Except for loading N x , load-concentration 
factors are somewhat less than those which exist in the structure with the six- 
member unit removed. The classical stress-concentration factors for an infi- 
nitely wide isotropic plate with a rigid inclusion, as listed in table I, are 
obtained from reference 9 with Poisson's ratio taken to be 1/3 to be consistent 
with the present studies. For extensional and shearing loadings, this classical 
stress-concentration factor is 1.50, and load- concentration factors for the mesh 
range from 1.28 to 1.37. Thus, classical stress-concentration factors over- 
estimate load-concentration factors for the mesh with the rigid insert, but not 
nearly as much as for the mesh with the hole. 


Models With Hexagonal Hole Patterns 

Since the load- concentration factors for a mesh with a six-member hexagonal 
hole were lower than stress-concentration factors in an isotropic plate with 
a circular hole, this study was extended to examine hexagonal hole patterns 
that encompassed more missing members, even though occurrence of such a pattern 
of damage is unlikely. To proceed with the study beyond the six-member hole, 
only one-quarter of the full model is used because of the symmetry of the 
model. 

Load-concentration factor as a function of r^/^, where r^ is the 
"radius" of the hexagonal hole or radius of the circle that circumscribes the 
hexagon, is shown in figure 10. Data are shown for arbitrarily chosen full 
mesh sizes of 8 x 8, 1 6 x 1 6, 24 x 24, 32 x 32, and 40 x 40. Data for the 8x8 
mesh are only determined for values of r^/^ equal to 1 and 2. To avoid plot- 
ting points very close to one another, data for all of the mesh sizes at values 
of r^/S, equal to 1 and 2 are not presented in figure 10. The solid line seg- 
ment in the figure corresponds to results for an infinite mesh obtained by 
linear extrapolation of the results obtained for the two largest meshes similar 
to the extrapolation done in figure 8. One thing that is not apparent from 
figure 10 is that except for loading N x , the relative location with respect to 
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the hole of the maximum loaded member may change as the hole size is increased. 
However, for a fixed hole size, the relative location with respect to the hole 
of the maximum loaded member remains unchanged as the mesh is increased. 

Table III presents the load- concentration factors for all the cases, but 
the individual load plots for each combination of mesh number and hexagonal 
hole size are not presented. An example of load distribution for a hexagonal 
hole with radius of 4 member lengths for a 32 x 32 mesh is presented in fig- 
ures 11(a), 11(b), and 11(c) for the loadings, N x , N X y, and Ny, respectively. 
The influence of the hole size for the loading N x on the 32 x 32 mesh is 
depicted in figures 12(a), 12(b), 12(c), and 11(a) for ratios of of 1, 

2, 3, and 4, respectively. The member loads increase as the hole extends over 
more rings of elements, and the effect of the hole diminishes for members far- 
ther from the hole. This is also true for the other loadings. 


CONCLUDING REMARKS 

A numerical study was conducted to examine the behavior of large space 
structures which have sustained damage or which have stiff equipment attached. 

A particular concept, the tetrahedral truss, was considered and attention was 
restricted to its planar faces under a series of basic loading conditions. When 
one member was removed near the center of the planar face of the tetrahedral 
truss, the maximum load in a nearby member increased by roughly 25 percent for 
extensional or compressive loadings and roughly 35 percent for shear loadings. 

When an entire six-member unit about a particular joint at the center was 
removed, maximum loads increased by roughly 30 percent to 45 percent for exten- 
sional loadings, the variation being dependent upon orientation, and roughly 
55 percent for shear loadings. Similarly, when a central six-member unit was 
replaced by a rigid insert (simulating stiff equipment), maximum member loads 
increased by roughly 35 percent. In general, the load-concentration effects 
investigated herein were found to be significantly less than the effects 
observed in classical stress-concentration problems for infinite plates with 
holes or rigid inclusions, although the effects tend to increase as a hexagonal 
hole extends over more rings of elements. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 22, 1979 
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TABLE II.- NONDIMENS IONAL MAXIMUM MEMBER LOADS FOR STRUCTURE 


CONTAINING SIX-MEMBER HEXAGONAL HOLE 


Loading N x 


Loading N xy 


Loading N y 


Loading, 


Ny corrected 


Mesh size 


P 0 ° 


Nx* 


4 x 4 
8x8 
16x16 
24 x 24 
28 x 28 
32 x 32 
40 x 40 


1 . 371 2 
1 .1 863 
1 .1297 
1 . 11 93 
1 . 11 70 
1 .11 56 
1 .1139 


p -60° 

N xy* 

p 60° 

Ny£ 

-2.0000 

1 .1547 

-1 .6641 

.9095 

-1.5649 

.8583 

-1.5453 

.8480 

-1 .541 1 

.8456 

-1.5383 

.8439 

-1.5350 

.841 8 


*60° 

Ny& 

1 .01 04 
.8971 
.8535 
.8442 
.8422 
.8408 
.8393 


TABLE III.- LOAD-CONCENTRATION FACTORS FOR STRUCTURE 


CONTAINING HEXAGONAL HOLES 



f " 

Load-concentration factors 

for 



Hole size, 

rh /A 






Mesh size 

Loading N x 

Loading N X y 

Loading 

N y 

Loading, Ny 
corrected 

00 

X 

00 

1 

1.370 

1.664 

1.529 


1.554 


2 

2.325 

2.631 

2.997 


2.989 

16 X 16 

1 

1.304 

1.565 

1.470 


1.478 


2 

1 .752 

1 .975 

2.282 


2.286 


3 

2.281 

2.777 

3.330 


3.342 


4 

3.078 

4.043 

4.931 


4.955 

24 x 24 

1 

1.292 

1.545 

1.457 


1.462 


2 

1 .658 

1 .846 

2.1 61 


2.1 62 


3 

1 .985 

2.358 

2.881 


2.886 


4 

2.370 

3.019 

3.711 


3.720 

32 x 32 

1 

1.288 

1.538 

1.453 


1.456 


2 

1 .626 

1 .800 

2.118 


2.119 


3 

1.890 

2.217 

2.733 


2.736 


4 

2.160 

2.694 

3.352 


3.357 

O 

X 

£» 

O 

1 

1.286 

1.535 

1.451 


1.454 


2 

1 .612 

1 .779 

2.098 


2.099 


3 

1.848 

2.152 

2.665 


2.668 


4 

2.069 

2.550 

3.1 94 


3.1 98 
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Figure 3.- Basic loading condi 
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(b) Loading Ny. 


Figure 5.- Continued 
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(c) Loading, Ny corrected 
Figure 5.- Concluded. 
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(b) Loading N y . 


Figure 6.- Concluded 
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(b) Loading N xy . 
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(c) Loadings, N y and N y corrected. 

Figure 8.- Maximum nondimensional member loads for structure with six-member 
hexagonal hole as function of mesh size, r^/H, = 1.0. 


24 






(c ) Loading Ny . 


to 

•vl 


Figure 9.- Concluded 



(a) Loading N x . 

Figure 10.- Load-ooncentration factor as function of hexagonal hole size. 
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Figure 1 0.- Continued 
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Figure 10.- Concluded 
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Figure 11*- Concluded 


(a) Hexagonal hole extending over one ring of elements. 

Figure 12.- Nondime ns ional member loads in 32 x 32 mesh structure loaded in x-direction with 
hexagonal holes of various sizes at center. Only one quadrant shown. 



I 



(JL> 

Ul 


(b) Hexagonal hole extending over two rings of elements 
Figure 12.- Continued. 
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